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Abstract 

C/j ' We study mappings obtained as s-periodic reductions of the lattice 

(-H Korteweg-De Vries equation. For small s € N 2 we establish upper bounds 

on the growth of the degree of the numerator of their iterates. These upper 
C"| ' bounds appear to be exact. Moreover, we conjecture that for any si,S2 

that are co-prime the growth is ~ (2s\S2)~ 1 n 2 , except when si + S2 = 4 
where the growth is linear ~ n. Also, we conjecture the degree of the n-th 
iterate in projective space to be ~ (si + S2)(2s\S2)~ 1 n 2 . 

in 

^ ; 1 Introduction 

(N 



Integrable mappings are characterised by low complexity [TJ 118) . This idea 

culminated in the notion of algebraic entropy, introduced by Viallet and col- 

f***^ ' laborators [4j [5] [7] . Low complexity means vanishing algebraic entropy which 

corresponds to polynomial growth of degrees of iterates of the mapping. A 
first proof of such a polynomial bound on the degrees was given in [3] . In [2 
it was proven that foliation by invariant curves implies zero algebraic entropy. 
Examples show that degree growth is a better indication of integrability than 
singularity confinement [7] [9] , cf . the discussion in |12) . Recently the notion 
has been extended to lattice equations [13 US] and used to find new integrable 
models [5]. 

In practise, one calculates the growth of degrees d n of the first n iterates 
of a mapping. Then one guesses the pattern by fitting the generating function 
d( x ) = Yl d n x n with a rational function p{x)/q{x) = g(x) and the algebraic 
entropy lim„_y 00 log(d„)/n is obtained as the logarithm of the inverse of the 
smallest zero of q(x), see [16] . We present an elementary method that enables 
one to derive upper bounds for the growth of degrees. Our formulas exactly 
produce all degrees that we have been able to calculate. 



2 Outline 

We will perform s-periodic reductions of the lattice Korteweg-De Vries equation 

(Ul, m - Uj+l,m+l)(«J+l,m - Ul,m+l) = OL. (1) 

This corresponds to studying solutions that satisfy the periodicity condition 
ui.m = ui+ sit m+s 2 - We choose si and S2 < Si to be co-prime natural numbers. 
Under this assumption the lattice equation reduces to a single ordinary difference 
equation (OAE) of order q := s± + S2 (or, a q-dimensional mapping). For 
background on periodic reductions we refer to [101 113] . There are q initial 
values, which we denote by X\,X2, ■ ■ ■ ,X q . The OAE, or the mapping, can be 
used to generate a solution a: n gz, which are rational functions in the initial 
values. 

One aim is to find a formula for the degree of the numerator (or denominator) 
of x n , as a function of n. We set x n — a n /b n , and derive a system of two OAEs 
for a n and b n , which are polynomials in the initial values. By choosing b n = 1 for 
n = 1, 2, . . . , q the degree (i.e., total degree in the variables x\ = a±, . . ■ , x q = a q ) 
of the numerator of x n is given by d% — d 9 n . Here d v n denotes the degree of a 
polynomial p n , and g n is the greatest common divisor g n = gcd(a„,6„). First 
we obtain a recursive formula for d® = d n + 1. Then we look at the growth of 
g n . After a number of iterates a miracle occurs: any divisor of b n will divide 
9n+ q (<? 7^ 4). This statement has been verified for a range of periodicities s, 
but seems to be difficult to prove in general. Next, we find a recurrence formula 
for the growth of the multiplicities of divisors: a divisor of g n divides g n +i 
with multiplicity ti, where t is an integer sequence satisfying a linear recurrence 
relation. We define a new set of polynomials c n — b n /f, where / is the product 
of all divisors of b i<n with the right multiplicities as given by the integer sequence 
t. Multiplying by / (which is a product c i<rl 's) and taking the degree on both 
sides of c n f — b„, we find that d c n + (d c * t) n = d b n where * denotes discrete 
convolution 

(d * t) n +i = ditn + d%t n -\ H V d n t\. (2) 

Using the recursive formulas for d b and t we find a recursive formula for d c , which 
can be solved to find polynomial growth of degree 2. Moreover, we obtain the 
coefficient of the leading term: (2siS2) _1 . 

We also consider the projective analogues of these mappings. We introduce 
homogeneous coordinates and derive a polynomial mapping in g-dimensional 
projective space. Here, the aim is to find a formula for degree of the components 
of this mapping. The strategy is very similar as the above. Once one has 
a divisor Cj of certain components of the mapping one can derive a recursive 
formula for the multiplicities at higher iterates of the mapping. At a certain 
point these multiplicities are (miraculously) higher than expected, after which 
the growth can be described recursively again. As before, a convolution formula 
provides us with a recurrence for the degrees of the divisors. In this case the 
degree of the n-th iterate is given by the sum 1 + d c n _ x + d c n _ 2 + • • • + d c n _ . 



This growth can also be described recursively and the leading term is found to 
be (si + s 2 )(2sis 2 )~ 1 n 2 . 

The case s — (3, 1) is exceptional. Here the growth is linear ~ n, and 
the mapping is linearisable. We provide its explicit solution in terms of an 
interesting sequence of polynomials, see section [3.31 and the appendix. 



3 Growth of degrees of rational mappings 

We first illustrate our approach by considering a low dimensional example, tak- 
ings = (2,l). 

3.1 A low dimensional example 

We take initial values Xi,X2,Xz on a staircase as in Figure Q] 

X\ x 2 



£l_£2 £3 £4 



Xl X2 X 3 X4, x 5 x e 



X3 Xi x 5 x 6 x 7 x s 



Figure 1: Staircase with (1, 2)-periodic initial values (xi,x 2 , £3), solved to the 
right. 



The x n are rational functions of x\, x 2 , X3, a which can be calculated recur- 
sively using 

x n = P(x n -i,x n - 2 ,x n -3), (3) 

where P solves equation (TTJ) for u/+i iJ7l , 

ct 

Ul+i, m = P(ut ym ,Ul+i >m +l,Ul,m+l) '■= Ul,m+1 H ■ (4) 

Ul, m — Uj+l,m+l 

We write x n — a n /b n . The recurrence © yields the following recurrences 
for a, b: 

a n = fln-3fn - ab n -ib n -2b n -z (5a) 

b n = b n ^ 3 w n (5b) 

where w n — a n - 2 b n -\ — dn-ib n -2- We choose b\ = b 2 = b 3 = 1, so that a n and 
b n are polynomials in the variables (initial values) a\,a% and a 3 . Their total 
degree will be denoted d£ an d 4> respectively. From (5) it follows that the 
degrees are at most 

< = maxKti + d a n _ 2 + d£_ 3 , <£_! + 4-2 + 4-1,4-1 + 4-2 + 4-3). 
4 = max(4_ 1 + d a n _ 2 + 4-3,4-1 + 4-2 + 4- 3 )- 



Given the initial degrees d° = d h n + 1 = 1 (n = 1, 2, 3) we find that 



< = 


= c 


-l+<-2 + C-3-l. 


4 = 


= <- 


-1+4-2 + 4-3 + 1 



(6) 

are upper bounds for the degrees of a n and 6„, and d„ = d h n + 1 (n e N). The 
sequence d b comprises sums of tribonaccci numbers, cf. [21 seq. A008937]. 
Certainly, these sequences grow exponentially. However, there will be a lot of 
cancelations in x n — a n /b n due to common factors of a n ,b n . We will prove that 
the degree of the greatest common divisor 

g n = gcd(a n ,b n ) 

is sufficiently large to ensure that d„ — d 9 n grows polynomially. 

Suppose that f tk divides g^ with k € {n— l,n — 2,n — 3}. Then from (5) it 
follows that /'" divides <?„, where 

t n = t n -l + in— 2 + tn-3- (7) 

We define an integer sequence t by t\ = t<2 = £3 — 2 = and the above recursion. 
Such numbers t are called tribonacci numbers, cf. [TJ] seq. A000073]. Thus 
have the following: 

! 2 \9n^! t3+ '\gn +l , I GN. 

By direct calculation, using Maple and the recurrences (5), we find that the 
polynomial w^ divides g n +3 (n > 3). This implies that w% divides g n +i- We 
can now write symbolically 

h = Ci = 1, i = 1,2,3, 
bi = c^ i = 4,5,6, 
&7 = c 7 c\, 

~l 2 2 td. t-i 

b s = csc^^csc^c^, 

b 9 = c 9 4clcl = c 9 c\*c%ci s c%c^c t ic* r 2 c t 8 1 , 



n-1 
Cn 



\{cr\ (8) 



which defines polynomials c„. Taking the degree of both sides of equation j8| 
we find d b n = d c n + (d c * t) n where * denotes discrete convolution, see @. From 
this we infer, using the recurrence for t ([7]), that 

d n ~ d-n — dit n -i + ■ ■ ■ + d c n _ i t i + d„_ 3 t3 

= dJ(tn-2 + tn-3 + tn-i) + ' ' ' + ^_ 4 (*3 + *2 + *l) + <£_ 8 *3 
= (d C * t) n _ x + {d C * t) n -2 + (d C * *)n-8 + 2c£_ 8 

— « n _i u n -l T "n-2 u n-2 ' "n-3 + u n-3 



which, using the recursion for d b ©, shows that 

d c n = d c n _ x + d c n _ 2 -d c n _ 3 + l. 

Together with d\ — d? 2 = d§ = this gives a sequence of quarter-squares, cf. 
QH seq. A033638], 

Note that the Cj< n 's in (|8]) are divisors of g n . Thus the quantity d„ — d^ is 
bounded from above by d b n + 1 — (d c * £)„ = d° + 1, which grows asymptotically 

~n 2 /4. 

3.2 More general periodic reductions 

Next, we consider the mapping obtained from s-periodic reduction taking s± and 
S2 to be coprime. Without loss of generality we may assume si > S2- Remember 
we denote s\ + S2 — q- Initial values X\,X2, ■ ■ • ,x q are given on a standard 
staircase [13] , see also [TO] in which a general theory of periodic reductions for 
equations not necessarily defined on a square has been developed. The initial 
values are updated by an recurrence of order q: 



(-71+1 



i\Xn — S2 t <En—si ) Xn—q) 5 



(9) 



cf. equation (j4}. For example, when s = (3, 2) we pose initial values as in Figure 
[2 These are updated by shifting over (2, 1), e.g. x§ >->• x 6 = P(x4,Xs,Xi). 

^i #3 x 5 



X2 Xi xe x$ 



X\ X3 X 5 X7 Xg Xn 



XA X 6 Xs XlO *12 Xu 

• o o o o o 

Figure 2: (3, 2)-periodic initial value problem updated in direction (2, 1). 



By setting x n — a n /b n we derive, for n > q, 



^n—q^n O^0 n _ Sl n _ S2 n _ 



®n On—q'uJn 



(10) 

(11) 



where 



On — si ^n- 



2 6 n _ Sl . We choose hi = 1, i = 1,2, ...,q so 



that a n and 6 n are polynomials in ai, a2, • • ■ , a q . As before, from initial degrees 
d° = d b n + 1 = 1 (n = 1, 2, . . . , tjr) we find that d° = rf* + 1 (n € N), and that 



d a = d a + d a +d a - 1 d = d + d +d +1 



are upper bounds for the degrees of a n and b n . If f tk divides g^ with k < n, 
then /*" divides g n , where 

In ^n— si i In—si i t n — q. l-L^J 

If initially t% = 0, i = 1, 2, . . . , g — 1, t g = 2 then 

Conjecture 1 TTie polynomial w^ divides g n + q (for n > q). 

It turns out this conjecture is more difficult to verify for s 2 << *i- We ver- 
ified the conjecture in the following ranges of values s 2 < si: s% = 1, . . . , 5 with 
s 2 <s 1 < 9s 2 , and s 1 = s 2 + l with s 2 = 6, 7, . . . , 25, 50, 100, 150, 200, 250, 1000. 

The conjecture would imply that w 1 ^ divides g n +i- Assuming it, we can 
define polynomials Ci by 



b n = c n Yl 



n-l 

t 
C 



1 1 



which yields d b n = d c n + (d c *t) n . Using the recurrences for t and d b we find 
d c = d c + d c - d c +1 

"n "n— si ' n-s 2 u n~q ~ ± - 

In the case s = (3, 2) the sequence [TH seq. A001399] 

0, 0, 0, 0, 0, 1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, ... , 

is given by 

JC 47 (-1)" C + C" n 1 1 a x» i 

gL c , = h — + n H n 2 , C = 1. 

" 72 8 9 2 12 ' s 

In general, the quantity d^ — d^ is bounded from above by d„ + 1, whose 
asymptotic growth is 

- (2sis 2 ) _1 n 2 . 

3.3 The exceptional case 

The case s = (3, 1) is an exceptional case. Here the growth is linear, which 
resembles the fact that that the mapping can be linearized. Introducing h = 
(xi — X3)(x 2 — xa), the mapping 

(X 

(x 1 ,X 2 ,X 3 ,X 4 ) l-> (X 2 ,X 3 ,X4,X! H ) 

x 4 - x 2 

reduces to h \— > a — h, which is an involutional Nevertheless, it is interesting to 
see what cancelations cause the growth to become linear. 



lr The function h is a 2-integral of the mapping. In |11| fc-symmetries are used to perform 
explicit dimensional reduction of mappings related to (si, l)-periodic reductions of lattice 
KdV. The dimension si + 1 is reduced to s± or s± — 2, when s± is even or odd, respectively. 



We set x n — a n /b n to find 

a n = a„_ 4 (a„_ 3 &„_ 1 - a n _i&„_ 3 ) - a&„„i&„_3&„_ 4 , (13) 

b n = b n -4,{a n -3b n -i - a n -ib n -a). (14) 

Taking initial values (01, 02, 03, 04) = (x — w, y + z, —w, z) and b\ = b 2 = &3 = 
64 = 1 we have found that, see the appendix, 

a„ = y u - 2 (a - sy)*»- 8 x*»-*c„, (15) 

b n = y 8 "+ 1 (a-xy) 8n x 8 "-\ (16) 

where to = ti = -s = 5i = and 

71 
*n+2=*n+l+*n+L4J> ( 17 ) 

ft 

Sn+2 = Sn+1 + S n + (-l) n LTj- 

Define r n = s„+4 — t n+ \. One can show that r„ = n(\ + (— l) n )/4, which 
is nonnegative. It follows that the a n /c n is a common divisor of a„ and 6„. 
Dividing out this factor we are left with denominator growth (n > 4) 

d b n ~ d a J c = r n -3 + 2r„_ 4 + r„_ 5 = n - 4. 

Note that in this case the common divisor of a n and b n consists of three different 
factors only, whereas for other values of s the number of common divisors grows 
linearly with n. Here, the multiplicity grows faster than what can be expected 
from the form of the recurrence. In other word a 'miracle' happens at every 
iterate: from (|17j) one can derive 

i n i 

tn+4 = tn+3 + t n +l + t„ + |_ — J , 

which should be compared to p^|) , taking si = 1, s 2 = 3, q = 4. 



4 Growth of degrees of projective mappings 

The entropy of a rational mapping has also been defined in terms of the growth 
of the degree of its equivalent in projective space J3J. Again we first consider 
the case s = (2, 1). 

4.1 A low dimensional example 

The 3-dimensional mapping is 

(xi,x 2 ,x 3 ) H> (x 2 ,x 3 ,xi H )■ 

X3 ~ X 2 



We set Xi — di/di, i = 1,2,3. If we denote the image by 6^/64, then the 
homogenised mapping is a i-> b: 



(18) 



Note that the first, second, and fourth component of the image share a common 
divisor. We are interested in the growth of the multiplicities of such a divisor. 
Suppose that c divides a\,a 2 and 04. From (|18[) it follows that c is a common 
divisor of 61, 63, 64. We continue the argument, 



a 2 


1— y 






/ a 2 (a 3 ~a 2 ) \ 
03(03 ~ o 2 ) 


03 
V a 4 J 








ai(a 3 - a 2 ) + ao^ 
\ a 4 (a 3 - a 2 ) / 



and 



c|(ai,a 3 ,a 4 ) =4> c|(6 2 , 63,64) 



c|(a 2 ,a 3 ,a 4 ) =>• c 2 |(&i, 6 2 , 64), c\b 3 



However, if we denote the common divisor of bi,b 2 , 04 by c, then miraculously 
c 3 divides all four components of the fourth iterate of a H> b. At the next 
iterates the multiplicities double. Denoting the multiplicity of c in the fourth 
component of the ith iterate by tj, we have t\ = t% = t 3 = 1, £4 = 3, and (at 
least) i n >4 = 2t n _i. We now introduce two sets of polynomials d, di as follows: 



a 2 

03 

\ o 4 ) 



m- 



H^ 



/ a 2 ci \ 

a 3 ci 

di 
\ a 4 ci / 

/ d 2 cfc3C4 

d 3 clc 2 C4 
d±c\c 2 c 3 
\ a 4 clc 2 c 3 C4, J 

( d„- 2 nr=i 3c * 
d«-i nr=i 2 

^ rUi C, + 

v 04 n l= i c, + 



\ 



/ a3Cic 2 \ / dic 2 c 3 
d 1 c 2 d 2 cic 3 

d 2 ci d 3 cic 2 

\ a 4 cic 2 J \ a±c\c 2 c 3 ) 

I d 3 c\clacz \ 

d 4 cfc2C 3 C5 

d 5 c\clc 3 Ci 
\ a 4 clclc 3 c 4 c 5 J 

\ 



H^ 



H^ 



h-> 



2 c *»+i 






(19) 



/ 



As an ordinary polynomial map the degree of the nth iterate is 

2 n = l + (d c *t) n+1 . 

Substracting 2™ = 2 + 2(d c *t) n from this equation, and using the recursion for 
t, we find 

d c n = d c n _ 1 + d c n _ 2 -d c n _ 3 + l. 

Projectively, the nth iterate (with n > 2) is 

(dn— 2 C n — lCn, d n —\C n — 2 C n , d n C n — 2 C n —i, a 4 C n — 2 C n —\C n ), 



after division by the common factor Yl7=i c i" +1 ■ We define 

n— 1 

Pn ■= a A 



n ^ 

i— max(l.n — 3) 



The projective degree is 



a n>3 



1 + d^_i + d° n _ 2 



*ro— 3' 



The recursion for d c yields c^ = d^_! + d p n _ 2 — d^_ 3 + 3. Together with initial 
values d? = 2 i_1 , i = 1,2,3 this gives the sequence [HI seq. A084684] 

1,2,4,8,13,20,28,38,49,62,..., 

which agrees with computations in projective space. The growth 

m 15 (-1)™ 3 3 2 

dP = h n + -n 2 

8 8 2 4 

is the same as for a mapping connected to the discrete Painleve I equation [4 



4.2 More general periodic reductions 

Now we consider the projective mapping that corresponds to s-periodic reduc- 
tion with si and s 2 coprime. We take si < s 2 , and q = s 1 + s 2 . It is convenient 
to take initial values Xq,x 2 , ... , X q -\. They are updated using the recurrence 
([5]), or equivalently, the g-dimensional mapping 

(jCo,Xl, ■ ■ -,Xq-l) !->• (»1, ■ • .,Xq-l,P(x ai ,X aa ,Xo))- 

Denoting the image of X{ = dija q by bi/b q we find a mapping a M> 6 in q- 
dimensional projective space, 



/ «o \ 

Ol 

a q -i 



( a\{x si -x S2 ) \ 

0-2{ x si — X S2 ) 

a (x Sl -x S2 ) +a 2 q 

\ a q\X Sl x s 2 ) J 



As in the case s = (2, 1) there is a common factor dividing all components but 
one. When s\ > 1 we have 



c|(o ,Oi,...,O g _2,Og) => 

When si > 2 we have 

J c 2 |(a ,ai, . . .,o g _ 3 ,o 9 _i,a g ), 
I c|a„_2 



c 2 \(b ,bi, . . . ,b q - S ,b q -i,b q ), 



c\b q - 2 



c 4 |(6 , bi, . . . , b g _4, b g -2, b g _i, b g ), 

9-3 



c 3 |fo. 



This doubling in most components continues until after si — 1 iterations we are 
lead to (if S2 > si + 1) 



c \{a ,ai,. 

2 al ~ 1 -l| 
c 0», 



■ ; o S2 _i, a S2 _|_i, . 




,b S2 -2,b s 



■,b q ), 



Then we have doubling again, until after S2 — 1 iterations where the growth is 
similar to the above. Doubling continues until ... 

Conjecture 2 The 'miracle' happens after q iterations where suddenly the mul- 
tiplicity is one higher than double the previous one. 

This we have only verified for a couple of small values of s±, S2- Conjecture [2] is 
harder to verify, using direct calculation, than conjecture Q] We will assume it 
in the sequel. We define integer sequences by ii = 1 and 



zt n 



n = si,s 2 , 
n = si + 82, 
otherwise. 



We now introduce two sets of polynomials a,di as follows: 



f a ° \ 


1 


Oi 






M> 


a q -2 




a q -i 




\a q J 


I 



/ a^ \ 



*2 Jl 



*1 

a 2 c{ 



dq-lCi 

dA 1 ' 1 



V agC* 1 



/ CL 2 C{ 

a-sc^c 1 ^ 



\ 



1 tx 



7 C 

aic^ c>2 

J *2 *1-1 

a 2 c^C2 
\ a q c 1 c 2 



H^ 



/ d lC \ q 



d 2 c 1 q c 2 q 






d q -lCi ■ ■ ■ C, 
t 



d q c x 



1 h 

9-1 °5 



tgl^ 






^ 



/ 



dsCi 






*„- 



9-1- 



2 °3 






V a q c x q 

( d n ^ q+1 nr=i c- ,i+i_ 7 c n-«+i \ 



. .. r *2 r *l-l 

C 9 C 9+l 



-tn-g+2 . 



-n-9+2 



in 



7 c n-i 



/ 



\ a g IL=i c t + 
As an ordinary polynomial map the degree of the nth iterate is 

2 n = l + (d c *t) n+1 . 



!->• 



H» 



(20) 



10 




Substracting 2 n = 2 + 2(d c * t) n from this equation, and using the recursion for 
t, we find 

1 < n < Si, 

si < n < s 2 , 
s 2 < n < q, 
q<n, 

or, d c n = d c n _ 8l + d c n _ S2 - d c n _ q + 1 for all n, taking d c n<1 = 0. 
Projectively, the last component of the (n — l)st iterate is 

n-l 

Pn ■= a q Yl C *' 

i— max(l,7i— q) 

which has degree 



We find 



c = i + E 

i— max(l,n- 


d\. 

-1) 


n, 

d p n _ Sl +n-l, 

dn- Sl + d l~s 2 +n-2, 

d p +d p - d p +a 


1 < n < si + 1, 

si + 1 < n < s 2 + 1 
s 2 + 1 < n < q, 
q < n. 



^>-{ / S '±/ 'j_„_0 c „x1/^„ ( 21 ) 



For example, in the case s = (2, 3) the sequence of degrees 

1, 2, 3, 5, 8, 12, 16, 22, 28, 35, 43, 52, 61, 72, 83, 95, 108, 122, 136, . . . 
is given by 

m 127 (-1)™ c,n-i +( -x-n 5 5 

d p = + - — ^ n + —n 2 , C = 1- 

72 8 9 6 12 ' s 

In general, the recursion (|2"Tj) yields asymptotic growth 

- {s 1 + s 2 )(2s 1 s 2 )- 1 n 2 . 

5 Conclusion 

In |17j Viallet discussed two approaches: the heuristic method, where no proofs 
are obtained, and serious singularity analysis, which is limited to 2-dimensional 
maps, or some exceptional higher dimensional cases. The question was raised 
how can we go further, in particular to high dimensions? One suggestion was 
given: the arithmetical approach, cf. 6 . In this paper we have presented a 
different approach and showed that it works for high dimensions, at least for 
(most) mappings obtained as reductions from an integrable lattice equation. 
The only condition on the dimension is that one has to be able to iterate the q- 
dimensional map q times to verify conjecture [1] or [2] The scope of this approach 
is left open for future research, e.g. to consider other reductions, other lattice 
equations, and non-integrable or almost-integrable maps. 



11 



A Solution of the (3,l)-map 

We prove that the recurrences (|13ll4p yields expressions (|15I16[) , with 

C2n+i = y(x - w){xy) n - 2 - P n , c 2n+2 = (z-y)(a- xy) 11 ^ 1 - yP n , (22) 
where 

n-l 

P n :=Y,T^ k {xy) h a n - 1 - k , 



fc=0 

with 

l k+l ~ ± k ~ ± k ' i ± Q - ± n~ - L - 

that is, [H seq. Al 12468] 



T'n+l rpn rpn-\-l rpn rpn -i 
7.11 — J- h -L i. , in — ^ „ — 1 , 



rpn.k _ V~V_i\ra-J| H ~*~ * ^ 

^-^ \ n — i 

Proof: Substituting (I15I16|) in (JT4J) yields 

C2n = (« -a;y)c 2 n-2 -y(xy) n ~ 2 , c 2n +i = xyc 2n -i - (a - xy) n ~ x . (23) 

Substituting (|22J) in (@3]) yields 

P n = (a- xy)P n -i + (xy)"- 1 , P n = («|/)^ n -i + (a - a^)" -1 

which can be verified using the definition of P and T. Substituting (|15I16|) in 
(fl3j) yields 



(xy)' 4 y(c 2i+ i - (a - xy) 1 2 ) = -c 2 i- 3 (c2 4 - (a - xy)c 2l - 2 ) 

and 

(a - xy) l ~ 4 (c 2l - a(xy) l ~ 3 y) = -c 2l - 4 (c 2l -i - xyc 2i - 3 ), 

which follow as a consequence of (f2"3"|) . □ 

Remark 1: The expressions for x n = a n /b n can be simplified as follows. Let 

(*i, £2)2:3, £4) = (a; - w,y + z,-w, z), x„> 4 = x„_ 4 H • (24) 

%n— 1 *£n— 3 

Then 

P P 

£2n+i = a: - w - a; ^ , x 2 n+2 = y + s - y (q _ gy)w _ 1 ■ (25) 

Remark 2: The recursion (IM1) can be solved explicitly as follows. The variable 
Vn = x n - x n+2 satisfies (y„ + y n -2)Vn-i = a and we find 

x\ - 2j=i V2i-i, n odd 
x 2 ~ l^iLi V2i, n even 
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We can solve y n in terms of z n — y n y n +i, 



„ _ ) Zn-2 Z„_ 4 Zl 9 " 

-2— --2—2 . . . =a.y 2 i n even 

(^ Z„_2 Z„_ 4 Z 2 yz ' 

and we have z n = a — z n -\ which implies 

Z\, n odd 

a — Z\ , n even ' 

where z\ — y\y<i — (xi — Xs)(x2 — £4) = xy. Backsubstituting yields 

(n-l)/2 




and we find (f2"5")) with 



P„ = g> - xy)'(«yr^ = ^^ . (26) 

* — ' a — Ixy 

Remark 3: From the latter expression (|26|) for the polynomials P„ it is easy 
to extract the following properties 

* They form a divisibility sequence: n\m =4» P n |P m . 

* The real parts of their zeros equal a/2: P n (xy) = ^ xy + xy = a. 
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